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Methanol masers in environments of three massive protostars 
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Abstract. We present the first EVN maps of 6.7 GHz methanol masers of three high-mass protostar candidates selected from the 
Torun unbiased survey of the Galactic plane. A variety of linear and arc like structures was detected. A number of maser clusters 
with projected sizes of 20-100 AU show monotonic velocity gradients. Some of them are roughly perpendicular to the major 
axes of these structures and can arise behind shock fronts. 



1. Introduction 

Observations of the 6668.519 MHz methanol maser transition, 
first detected by Menten ( 1991 1, appear to be powerful tools 
to identify the massive early type stars still embedded in their 
parental dense molecular clouds. The high brightness of this 
line enable us to investigate structures at milliarcsecond (mas) 
scales (a few hundreds of AU at the distances of a few kpc). 

Phillips et al. (1998! analyzed 45 methanol maser sources 
in star-forming regions and divided them into five groups on 
the basis of their morphology: linear (curved), elongated, pair, 
complex and simple. The linear masers were outstanding in 
their survey and showed a monotonic or near-monotonic ve- 
locity gradient along the source major axis that is consistent 
with a model of masers embedded in a rotating disk (radius 
of a few thousands of AU) seen edge-on or nearly edge-on 
around a high mass (up to 120 M Q ) protostar or OB star (Norris 
et al. ITffJgl Minier et al. EOOOV However, Walsh et al. (1998 1 
found only 12 sources showing velocity gradients in the sam- 
ple of 97 methanol sites studied. They proposed a scenario in 
which the methanol masers appear before the UC H II phase 
around the protostar, associated with embedded non-ionizing 
stars. The dense knots of gas are compressed and accelerated 
by the passage of the shock and local conditions cause different 
geometries of maser spots. Dodson et al. (2004i, based on the 
VLBl data, proposed that the 6.7 GHz methanol masers arise 
behind low-speed (<10kms~') planar shock. A shock propa- 
gating nearly perpendicular to the line-of-sight produces a lin- 
ear spatial distribution of maser components. Interaction of the 
shock with density perturbations in the star-forming region dis- 
rupts the linearity of maser structures. 

The VLBl results presented below are first in a series for 
a large sample of methanol masers detected in the unbiased 
survey of the Galactic plane (Szymczak et al. 2002 1. Our aims 
are twofold; to investigate the mas scale of 6.7 GHz methanol 
maser structures in order to test the above mentioned hypothe- 
ses and to search for relationships of methanol emission with 
other tracers of star-forming activity. We improved the posi- 
tions of newly detected 30 sources in the Torun survey using the 
single baseline of MERLIN (Mark II and Cambridge antennas). 



Table 1. Coordinates (J2000) of the target sources. 
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G33.64-0.21 


18 53 32.551 


00 32 06.525 


0.3 


4 


G35.79-0.17 


18 57 16.911 


02 27 52.900 


0.6 


3 


G36.11+0.55 


18 55 16.814 


03 05 03.720 


0.2 


1.4 



The absolute positions of the strongest maser features were de- 
termined with accuracies of 0.1-078 in RA and 1-4" in Dec 
depending on the source intensity and declination (Bartkiewicz 
et al., in prep.). 



2. Observations and data reduction 

The observations of methanol maser emission at 
6668.519 MHz of the three star-forming regions G33.64-0.21, 
G35.79-0.17 and G36.11+0.55 1 were carried out with the 
European VLBl Network (EVN) on 2003 June 08. Useful data 
were obtained from four antennas (Cambridge, Effelsberg, 
Jodrell Bank and Onsala). In the single dish survey these 
sources showed complex methanol spectra with three of more 
clearly visible features of flux densities higher than 10 Jy. All 
of them are OH emitters (Szymczak & Gerard 2004 1. The 
coordinates of the target sources and their errors as measured 
with the Mark II - Cambridge baseline of MERLIN are given 
in Table 1 . 

The targets were observed for a total of 12 hr, together with 
observations of continuum source 3C345 for the purpose of 
bandpass, delay and rate calibration. The source J1907+0127 
(0.2 Jy at 6.7 GHz) was used as a phase calibrator for all three 
targets. The cycle time between each target and the phase- 
calibrator was 5.5 min + 3.5 min. We used a spectral bandwidth 
of 2 MHz, resulting in the velocity coverage of 1 00 km s~ 1 , di- 
vided into 1024 channels at correlation to give a velocity res- 
olution of 0.09 km s~'. The bandwidth was centred at the local 
standard of rest (LSR) velocity of 65 km s . Observations were 
made in left- and right-hand circular polarization, but since 
the methanol emission is not strongly circularly polarized both 
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polarizations were averaged in order to improve the signal-to- 
noise ratio. 

The data were correlated on the EVN Mk IV Data Processor 
operated by JIVE. We carried out the data reduction with 
standard procedures for spectra line observations using the 
Astronomical Image Processing System (AIPS). To make final 
maps we used 1 mas pixel separations and an elliptical restor- 
ing beam 14x6 mas with a PA of -1° and applied the uniform 
weighting. The rms noise level was of 7-10mJy beam -1 in 
emission-free Stokes / maps. 

We created fringe rate maps of the brightest channels of the 
targets but still we failed to find the absolute position of our 
sources. The target sources were near zero declination (from 
+(K5 to +3°.l). The phase-calibrator J1907+0127, the nearest 
one from the VLBA Calibrator Survey, was about 3° apart from 
all three targets. Furthermore, due to use of only four EVN tele- 
scopes listed above the uv-plane coverage was poor for N-S 
baselines. It is likely that these factors together with too long 
phase-referencing cycle time precluded a proper phase calibra- 
tion. However, the positions of the strongest components of 
each target are known with accuracies better than (X'6 in RA 
and 4" in Dec (Table 1). 



3. Results and Discussion 

3.1. G33.64-0.21 

Fig. 1 shows the spectrum and the overall distribution of maser 
emission in G33.64-0.21. The emission is seen in the veloc- 
ity range from 58.5 to 63.3 km s~' and 94 maser components 
brighter than 40 mJy b -1 (5cr) are detected. These form 12 clus- 
ters distributed along an arc like structure of overall size of 
about 150 mas that corresponds to 585 AU for the assumed 
near kinematic distance of 3.9 kpc. We do not note any regular- 
ity in the velocity of the whole structure that can be expected 
for a rotating disk. However, monotonic velocity gradients are 
found within six clusters named as A, C, F, G, J and K (Fig. 1). 
Angular sizes of individual clusters range from 5 to 18 mas i.e. 
their projected linear sizes range from 20 to 70 AU. We notice 
that for most clusters the velocity gradients are roughly per- 
pendicular to the arc. We state, the methanol maser emission in 
G33.64-0.21 does not show large-scale velocity gradient but 
about half of the maser clusters exhibit internal velocity gradi- 
ents some of which are roughly perpendicular to the arc. This 
can suggest that the maser structure forms in a shock front pos- 
tulated by Dodson et al. (2004). 

The inset in Fig. 1 shows that amplitude of feature L of 
the EVN spectrum is comparable (within 30%) with that of the 
single dish spectrum. For the rest of features the amplitudes 
observed with the EVN are a factor of 3-4 lower. This implies, 
even though variability and/or errors in flux density cannot be 
ruled out, that a large fraction of maser flux was missed in the 
interferometric observations. Therefore, G33.64-0.21 appears 
as a good candidate to map a low intensity and diffuse maser 
emission. 



3.2. G35.79-0.17 

The spectrum of G35.79-0.17 is composed of four features 
within the velocity range from 59.9 to 62.7 km s~' (Fig. 2). 
The shape of the spectrum is the same as observed in 2000 
(Szymczak et al. 2002) but amplitude is about a factor of 2.5 
lower. Variability of the source and calibration errors cannot 
fully account for such a difference. We suggest that there is dif- 
fuse emission missed during our EVN observations. We found 
the emission in 33 spectral channels and its distribution ap- 
pears as a 9.5 mas linear structure (Fig. 2). Analysis of the 
position-velocity diagram revealed a linear velocity gradient 
along NE (red-shifted) - SW (blue-shifted) direction. Such a 
characteristic structure can arise from masers from an edge-on 
rotating disk or accelerating outflow. In a case of an edge-on 
rotating Keplerian disk we can estimate the lower limit of the 
mass of a central star (e.g. Phillips et al.[F?98). Assuming the 
far and near kinematic distances of 4.6 kpc and 10.3 kpc we 
infer that the enclosed mass is 0.4 and 0.9 M , respectively. 
These values are certainly underestimated; we have assumed 
that there is no inclination and that the disk size equals the 
extension of the maser emission. Minier et al. (2000 1 found 
similar structures in W51 and G 29.95-0.02. They also de- 
rived sub-solar masses and concluded from a larger sample 
that they had detected only a small fraction of the disk, which 
typically extends over 1000 AU. Indeed, the maser structure in 
G35.79-0.17 is similar in a size with individual clusters ob- 
served in G33. 64-0.21 (e.g. the cluster A). We cannot exclude 
that the observed structure arises due to an outflow from a mas- 
sive protostar signposted by infrared source IRAS 1 8547+0223. 
Proper motion studies would help to solve the problem of ori- 
gin of the linear structure of methanol maser in G35.79-0.17. 

3.3. G36.11+0.55 

Fig. 3 shows the spectrum and the distribution of the 6.7 GHz 
methanol maser emission in G36.ll +0.55. We found two ac- 
tive regions, north-western and south-eastern separated by 1'.' 1 . 
The maser components in both regions are elongated along a 
position angle of ~45°. These structures are nearly perpendic- 
ular to the major axis. The NW region is composed of eight 
clusters, labelled from A to H, within the LSR velocity range 
from 81.1 to 84.5 km s~'. The SE region is composed of six 
maser clusters, labelled from I to N, of the LSR velocity from 
73.0 to 76.2 km s~'. 4 out of 14 clusters exhibit monotonic ve- 
locity gradients; two of them (NW region) are perpendicular 
to the major axis. For the assumed near kinematic distance of 
5.3 kpc the two clumps NW and SE are separated by 5840 AU. 

The methanol emission of G36. 11+0.55 lies 5'.'5 SE 
from the centre of a giant molecular clump Mol77 
(IRAS18527+0301) with embedded high-mass protostars of 
type B2.5-08.5 (Brand et al. 1200 1> . Since the absolute po- 
sition of methanol masers is known within 0V2 and 1'.'4 (RA 
and Dec) this positional offset may be significant. The veloc- 
ities of thermal emission of CS, CO and HCO + range from 
71.9 to 82.1 kms -1 with the peaks about 75.5 km s _1 and match 
well the velocity range of methanol emission. We note that 
the overall structure delineated by methanol masers is nearly 
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Fig. 1. The EVN 6.7 GHz methanol maser distribution in G33.64-0.21. The coordinates are relative to the strongest maser 
component (Table 1). Each triangle corresponds to a single component found in the channel maps. Maser clusters A, C, F, G, J 
and K show clear velocity gradients traced by arrows (pointed to the blue-shifted velocities). The spectrum (inset) is composed 
of emission profiles of maser clusters labelled from A to L. Thin line shows the single dish spectrum taken in 2000 and scaled 
with the right hand ordinate. 
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Fig. 2. The EVN 6.7 GHz methanol spectrum (left) and the distribution of maser components (right) in G35.79-0.17. The co- 
ordinates are relative to the brightest maser component (Table 1). Each circle corresponds to a single component found in the 
channel maps. The sizes of circles are proportional to the logarithm of the brightness. The numbers indicate the LSR velocities 
(kms -1 ) of the components. 
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Fig. 3. The overall distribution (top right) and spectrum (bottom right) of the 6.7 GHz methanol maser emission in G36. 1 1 +0.55. 
The coordinates are relative to the strongest maser component (Table 1). Enlarged areas of the NW and SE emission are shown 
on the left side. Each triangle corresponds to a maser component found in the channel maps. Arrows mark velocity gradients 
similarly as in Fig. 1. The spectrum shown by thin line is the sum of spectra of individual maser clusters (thick line). 



perpendicular to the axis of an outflow traced by the blue- 
(71.0-73.8kms 1 ) and red-shifted (79.5-83.2kms" 1 ) emis- 
sion of the 13 CO line as reported by Brand et al. (2001 1. 
Therefore it is very likely that the 6.7 GHz maser emission 
mapped with the EVN are physically related to a high mass 
protostar. 

4. Conclusions 

The 6.7 GHz methanol maser emission towards three star- 
forming regions was imaged with milliarsecond resolution. 
Linear and arc like structures were detected. There were con- 
siderable numbers of maser clusters with internal velocity gra- 
dients roughly perpendicular to the major axis. These geomet- 
rical structures can arise in outflows and behind shock fronts. 
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